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Abstract 



The mathematical pendulum is the simplest system having all the basic 
properties inherent in dynamic stochastic systems. In the present paper we 
investigate the mathematical pendulum with the aim to reveal the properties 
of a quantum analogue of dynamic stochasticity or, in other words, to obtain 
the basic properties of quantum chaos. 

It is shown that a periodic perturbation of the quantum pendulum (sim- 
ilarly to the classical one) in the neighbourhood of the separatrix can bring 
about irreversible phenomena. As a result of recurrent passages between de- 
generate states, the system gets self-chaotized and passes from the pure state 
to the mixed one. Chaotization involves the states, the branch points of whose 
levels participate in a slow " drift" of the system along the Mathieu character- 
istics this "drift" being caused by a slowly changing variable field. Recurrent 
relations are obtained for populations of levels participating in the irreversible 
evolution process. It is shown that the entropy of the system first grows and, 

after reaching the equilibrium state, acquires a constant value. 
PACS number: 87.22.Jb 
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I. INTRODUCTION. FORMULATION OF THE PROBLEM 



Dynamic stochasticity is directly connected with the assumption that classical equations 
of motion may contain nonlinearities which arise when the (exponential) repulsion of phase 
trajectories occurs at a sufficiently quick rate. In the case of quantum consideration, the 
dynamics of a system is described by a wave function that obeys a linear equation, while the 
notion of a trajectory is not used at all. Hence, at first sight it seems problematic of find 
out the quantum properties of systems whose classical consideration reveals their dynamic 
stochasticity. A quantum analogue of classical stochastic motion is usually called quantum 
chaos. 

On the other hand, it is of practical interest to investigate parametrically dependent 
Hamiltinians H(Q, P, I), where (Q, P) is the set of canonical coordinates, / is the parameter 
describing how the system is related to the external field. The interest in such systems is 
explained by their use in the study of quantum points and other problems of mesoscopic 
physics [1]. 

In most of the papers that deal with parametrically dependent systems, their authors 
consider the following situation. For / = 0, the Hamiltonian H(Q, P, 0) is exactly integrable. 
As / increases, the Hamiltonian H(Q, P, /) becomes nonintegrable and, for a certan value of 
l , solutions of the classical equations corresponding to H(Q, P, l ) become chaotic. In the 
case of quantum consideration, eigenvalues E n (l ) and eigenf unctions ip n (lo) are found in the 
above-mentioned area of parameter values by using the method of numerical diagonalization. 
In that case, we show interest in the dependence of the parametrical kernel P(n/m) = \ < 
4>n{lo + ^OlV'm(^o) > | 2 on a parameter displacement 51 <C I. The value P(n/m), averaged 
statistically over states n, P(n/m) = P(n/n + r) = P(r) can be interpreted as the local 
density of states. The introduction of P{r) means that we pass from the quantum-mechanical 
description to the quantum statistical description [2 - 5] carried out by an intuitive reasoning. 

In the problems considered in the above-listed papers, the Hamiltonian H(Q, P, I) dis- 
plays chaos for both parameter values / = Iq and I = l + 51. 
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As different from these papers, in the present paper we investigate the situation, in 
which the Hamiltonian H(Q,P, I) is integrable and becomes nonintegrable after adding a 
strictly periodic perturbation 5l(t). As the basic Hamiltonian we take the Hamiltonian of 
the mathematical pendulum (universal Hamiltonian). 

As is known, the Schrodinger quantum-mechanical equation for the universal Hamilto- 
nian is written in the form of the Mathieu equation. The Mathieu-Schrodinger equation 
for an atom, which is under the action of optical pumping in the area of large quantum 
numbers, was for first time obtained by G. Zaslavsky and G. Berman [6]. These authors 
also performed analysis of quasiclassical states of the Mathieu-Schrodinger equation [6]. 

The main objective of the present paper is to investigate the behavior of the quantum 
mathematical pendulum in the area of dynamic stochasticity parameters. As is known [7], 
this area, called the stochastic layer, lies in the neighborhood of the separatrix of the classical 
pendulum. 

We show here that with the appearance of quantum chaos the pure state passes to the 
mixed one. In other words, the reversible quantum process transforms to the irreversible 
process of quantum chaos which can be described by a kinetic equation. The common feature 
of classical dynamic chaos and quantum chaos is, as will be shown below, the irreversibility 
of their states. 

II. A PARAMETRIC ALLY DEPENDENT HAMILTONIAN 

After writing the stationary Schrodinger equation 



Hi{> n = 



n 



(1) 



for the universal Hamiltonian of the atom+pumping system 



H = 



dip 2 



+ V, 



V = I COs2<y9, 
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we come to the equation coinciding with the Mathieu equation [8, 9] 



dip 2 



+ (E n - V(l, (p))i/j n = 0, (2) 



where E n — > jp^ are the introduced dimensionless values, / is the dimensionless amplitude 



of pumping, to' = ^jp- is the derivative of nonlinear oscillation frequency with respect 
to the action I [9]. 

The Mathieu-Schrodinger equation is characterized by a specific dependence of the spec- 
trum of eigenvalues E n {l) and eigenfunctions ip n (f, on the parameter / (see Fig.l). On the 
plane (E,l) with the spectral characteristics (so-called Mathieu characteristics [9]) of the 
problem, this specific feature manifests itself in the alternation of areas of degenerate (G T ) 
and nondegenerate (G) states (see [10], Figs. 3, 4). The boundaries between these areas pass 
through the branch points of energy therms E n (l). 

Degenerate and nondegenerate states of the quantum mathematical pendulum were es- 
tablished by studying the symmetry properties of the Mathieu-Schrodinger equation. In 
[10], by using the symmetry properties of the Mathieu-Schrodinger equation and applying 
the group theory methods, the eigenvalues for each of the areas G + , G_, G were found: 

G- -> = -y (ce 2n+ i(y?) ± ise 2n+1 (ip)), (3GL) 

G -> ce 2n (y?); ce 2 „ + i(y?); se 2n (y?); se 2n+1 ((p), (3G) (3) 

G+ -> Ctni'P) = ^(ce 2n (y?) ± ise 2n+1 (<f)), (3G+) 

C 2 n+i(v?) = ^|(ce 2n+ i(v?) ±ise 2n+2 (y?)). 

Here ce n (ip) and se n (<^) denote the Mathieu functions [9]. 

The wave functions (3G±) and (3G) form the bases of irreducible representations of 
the respective groups. Each of the four functions (3G) forms a one-dimensional irreducible 
representation of the Klein group V, while the functions ?p 2n+ i ((f) , 4> 2n ((p) from (3G_) and 



itniv) Atn+iiv) from (3G+) form the two-dimensional irreducible representations of two 
invariant subgroups of the group V [10]. 

Let us assume that the pumping amplitude is modulated by a slowly changing elec- 
tromagnetic field. The influence of modulation can be taken into account by making a 
replacement in the Mathieu-Schrodinger equation 

l(t) -> l + Al cos id, (4) 

where Al is the modulation amplitude expressed in dimensionless units, v is the modulation 
frequency. 

We assume that a gradual change of l(t) may involve some N branch points on the left 
and on the right side of the separatrix (Fig.l): 

Al>\l n + -l% n= 1,2,... N. (5) 

After making replacement (4) in the universal Hamiltonian, we obtain 

H = H + H'(t), 



dip 2 



Ho = — + l oCOs2ip, (6) 



H'(t) = Al cos 2cp cos vt. (6') 

Simple calculations show that the matrix elements of perturbation H'(t) with respect to 
the wave functions (3G) of the nondegenerate area G are equal to zero 

{ce n \H'(t)\se n ) ~ Al J ce n ((p) cos 2(pse n ((p)d(p — 0, (7) 

o 

where n is any integer number. Therefore perturbation (6') cannot bring about passages 
between nondegenerate levels. 

The interaction H'(t), not producing passages between levels, should be inserted in the 
unperturbed part of the Hamiltonian. The Hamiltonian obtained in this manner can be 
considered as slowly depending on time. 



Thus, in the nondegenerate area G the Hamiltonian can be written in the form 

l(t) = l + Al cos id. (8) 

As has been mentioned above, to different areas on the plane (E, I) we can assign different 
eigenfunctions (3G_), (3G), (3G+). Because of the modulation of the parameter l(t) the 
system passes from one area to another, getting over the branch points. 

III. PASSAGE FROM THE QUANTUM-MECHANICAL DESCRIPTION TO THE 
KINETIC DESCRIPTION. IRREVERSIBLE PHENOMENA 

As different from the nondegenerate states area G, in the areas of degenerate states G_ 
and G + , the nondiagonal matrix elements of perturbation H'(t) (6') are not equal to zero. 
For example, if we take the matrix elements with respect to the wave functions iptn+i ( see 
(3Cr_)), then for the left degenerate area GL it can be shown that 

H' + _ = H'_ + =< ^Vil^'WlVWi >~ AJ J ^ n+ i^: + i cos 2^ ^ 0. (9) 

o 

Note that the value H' + _ has order equal to the pumping modulation depth Al. 

Analogously to (9), we can write an expression for even 2n states as well. 

An explicit dependence of H'(t) on time given by the factor cos vt is assumed to be slow 
as compared with the period of passages between degenerate states that are produced by the 
nondiagonal matrix elements H' + _. Therefore below the perturbation H' + _ will be assumed 
to be the time-independent perturbation that can bring about passages between degenerate 
states. 

Thus, in a degenerate area the system may be in the time- dependent superpositional 

state 

iMO = c+(t)4>+ n + c-(t)^ n . (fo) 
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The probability amplitudes C^(t) are defined by means of the fundamental quantum- 
mechanical equation expressing the casuality principle [11]. We write such equations for 
a pair of doubly degenerate states: 

dC + 

-ih-jje- = (E on + H' ++ )Ct + H' + _C-, 
An- 

-ih-j± = H' + _C: + (E on + H'__)C~. (11) 

Let us solve system (11). In our case it can be assumed that H' ++ = H'_ and H' + _ = H'_ + . 

Let us investigate changes that occurred in the state of the system during time AT while the 
system was in the area G_, assuming that AT is a part of the modulation period T, AT < T. 
For arbitrary initial values system (11) has the solution 

C+(t) = ei Et (C + cos(^-t) + iC_sm(^-t)), 
C~(t) = ei Et (C_cos(^t) +i C+ sm(^t)), (12) 

where we redenote E a + H' ++ — > E, H' + _ — > H' . A slow dependence of the interaction H'(t) 
(6) on time can be taken into account in (12) if we use the replacement H' — > H' cos vt. 

Let motion begin from the state ip^n °f the degenerate area. Then as the initial conditions 
we take 

C-(0) = 1, C+(0) = 0. (13) 
Having substituted (13) into (12), for the amplitudes C^(t) we obtain 

Cn(t) = iexp(-Et) smut, 

I v 

C~(t) = exp(^Et) cos cut, (14) 

where the frequency of passages between degenerate states, r is the passage 

time. 

Note that the parameter u, which is connected with the modulation depth A/, has 
(like any other parameter) a certain small error Sou, which during the time of one passage 



t ~ 27r /uu, leads to an insignificant correction in the phase 2ti(5u/uj). But during the time 
t ~ AT, there occur a great number of oscillations ( phase incursion takes place) and, in 
the case AT ^> r, a small error Sou brings to the uncertainty of the phase ~ ATSou which 
may have order 2tt. Then we say that the phase is self-chaotized. 

Let us introduce the densitity matrix averaged over a small dispersion 5u: 



1 W+(t) iF n (t) ^ 
v -iF*{t) W-{t) } 



(15) 



where W^(t) = |C±(t)| 2 , F n (t) = C+(t)C~*(t). The overline denotes the averaging over a 
small dispersion 5u 



Ul+Su! 



A ^' t) = 2^j J A{ - X ^ dx ^ 

ui—Slo 



To solve (16) we can write that 



I. 



W+(t) = sin 2 cut, W~(t) = cos 2 cut, F n (t) = -sm2cut. (17) 

After a simple integration of the averaging (16), for the matrix element (17) we obtain 

^) = ^(1 T /(2M cos2wf), 

F n (t) = F* n {t) = ^f(25out) sin2cut, (18) 
„. . , sin25ujt 

At small values of time t <CT (r = 2-n/5u), insufficient for self-chaotization (f(25cut) ~ 
1), we obtain 

W+(t < r) = sin 2 out, W~(t < T ) = cos 2 out, F n (t < r) = - sin 2cjt. 

Z 

Comparing these values with the initial values (17) of the density matrix elements, we 
see that the averaging procedure (16), as expected, does not affect them. Thus, for small 
times we have 
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P+"(f «r) = 



sin ut | sin 2oot 
sin 2u;t cos 2 ut 



(19) 



One can easily verify that matrix (19) satisfies the condition p 2 n (t «r) = p n (t <C r), 
which is a necessary and sufficient condition for the density matrix of the pure state. 

For times even smaller than t <C r <C T, when passages between degenerate states 
practically fail to occur, by taking the limit out <C 1 in (19), we obtain the following relation 
for the density matrix: 



(t = 0) = p£-(* « = 



(20) 



This relation corresponds to the initial relation (13) when the system is in the eigenstate 
ip2n- Let us now investigate the behavior of the system at times t >T when the system gets 
self-chaotized. 

On relatively large time intervals t > r, in which the self-chaotization of phases takes 
place, for the matrix elements we should use general expressions (18). The substitution of 
these expressions for the matrix elements (18) into the density matrix (15) gives 



Ptr(t) = \ 



\ 



(21) 



1 - f(28ut) cos 2ut f(26ut) sin 2ut 
y -if(26ut) sin 2ut 1 + f(26ut) cos 2ut 

Hence, for times t > T during which the phases get completely chaotized, after passing to 
the limit 5uot 3> 1 in (21), we obtain 

f I- 0(e) iO(e) ^ 
v -iO{e) I + 0(e) J ' 



P+-(*»r) = | 



(22) 



where 0(e) is an infinitesimal value of order e = -^tjl' 



The state described by the density matrix (22) is a mixture of two quantum states tp^n an d 
ip2n with equal weights. The comparison of the corresponding matrix elements of matrices 
(22) and (21) shows that they differ in the terms that play the role of quickly changing 



fluctuations. When the limit is t ^> r, fluctuations decrease as ~ 

9 



(see Fig. 2 and 3). 



Thus the system, which at the time moment t — was in the pure state with the wave 
function ip 2n (20), gets self-chaotized with a lapse of time t ^> r and passes to the mixed 
state (22). In other words, at the initial moment the system had a certain definite "order" 
expressed in the form of the density matrix p + ~(0) (20). With a lapse of time the system 
got self-chaotized and the fluctuation terms appeared in the density matrix (21). For large 
times t ^> t a new "order" looking like a macroscopic order is formed, which is defined by 
matrix (22). 

After a halfperiod the system passes to the area of nondegenerate states G (20). In 
passing through the branch point, there arise nonzero probabilities for passages both to the 
state ce 2n and to the state se 2n . Both states ip 2n an d ^ 2n wm contribute to the probability 
that the system will pass to either of the states ce 2n and se 2n . For the total probability of 
passage to the states ce 2n and se 2n we obtain respectively 

1 1 2w 

PiPtnit > r ) Ce 2n) = d- / ^tni.V)ce 2n {.^)d^\ 2 + 

Z 71 J 

1 1 ^ 

+-\- J ^Mce 2n (^\ 2 2 2 2 - 2 . 

o 

2tt 

11/ 

P (pin(t > r ) -»■ S e 2n) = gl-y ^ + n(v)^2n(^)^|' 



2tt 

1 1 /• 2 11 11 1 



1 1 2lT 11111 



o 



Thus, in the nondegenerate area the mixed state is formed, which is defined by the density 
matrix 



A (t ~ ? » t) = i 



(24) 



where i and /c number two levels that correspond to the states ce 2n and se 2n . 

As follows from (24), at this evolution stage of the system, the populations of two nonde- 
generate levels get equalized. It should be noted that though the direct passage (7) between 
the nondegenerate levels is not prohibited, perturbation (6') essentially influences "indirect" 
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passages. Under "indirect" passages we undestand a sequence of events consisting a passage 
G — > G_ through the branch point, a set of passages between degenerate states in the area 
G-, and the reverse passage through the branch point G_ — > G. The "indirect" passages 
ocurring during the modulation halfperiod T/2 result in the equalization (saturation) of two 
nondegenerate levels. 

As to the nondegenerate area, the role of perturbation H'(t) in it reduces to the dis- 
placement of the system from the left branch point to the right one. 

It is easy to verify that after states (24) pass to the states of the degenerate area G + , we 
obtain the mixed state which involves four states ifni'-P) an d C^+ilv 9 ) ( see Fig-1)- 

Let us now calculate the probability of four passages from the mixed state pf n (22) to 
the states Ctni'P) and (t-M- 



1,1 

2 



2tt 



- / (ce 2 „(v?) + se 2n {^))if n {^)d^ 
ix J 



4' 



II 2 / 1 



(25) 



As a result of these passages, in the area G + we obtain the mixed state described by the 
four- dimensional density matrix 



P2„; 2 n+i(^ ~ T > r) 



V 



(26) 



/ 



^ 1 ^ 
10 
10 
1 

where the indices of the density matrix (26) show that the respective matrix elements are 
taken with respect to the wave functions ^tnif) an d C^i+ilv 9 ) °f degenerate states of the 
area G + . 

It is easy to foresee a further evolution course of the system. At each passage through 
the branch point, the probability that an energy level will get populated is equally divided 
between branched states. We can see the following regularity of the evolution of populations 
for the next time periods. 
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After odd halfperiods, the population of any n— th nondegenerate level is defined as an 
arithmetic mean of its population and the population of the nearest upper level, while after 
even halfperiods - as an arithmetic mean of its population and the nearest lower level. This 
population evolution rule can be represented both in the form of Table 1 and in the form of 
recurrent relations 

P[n, 2k] = P[n + 1, 2k] = \{P[n, 2k - 1] + P[n + 1,2k - 1]), 

P[n + l,2k + l] = P[n + 2, 2k + 1] = \{P[n + 1, 2k] + P[n + 2, 2k]), (27) 

where P[n, k] is the population value of the n— th level after time k^, where k is an integer 
number. The creeping of populations among nondegenerate levels is illustrated in Fig.4. 

The results of numerical calculations by means of formulas (27) are given in Fig. 5 and 
Fig. 6. Fig.5 shows the distribution of populations of levels P(n) after a long time t > T 
when the population creeping occurs among levels, the number of which is not restricted 
by (5). Let us assume that at the initial time moment t — 0, only one n — th level is 
populated with probability P(n ) = 1. According to the recurrent relations (27), with a 
lapse of each period T "indirect" passages will result in the redistribution of populations 
among the neighboring levels so that, after a lapse of time t = kr 3> T, populations of the 
extreme levels will decrease according to the law [10] 

P(n ±fc)~l 

If the number iV of levels defined by condition (5) is finite, then, after a lapse of a long 
time, passages will result in a stationary state in which all N levels are populated with the 
same probability equal to 1/N (see Fig. 6). The distribution obtained by us is analogous to 
the distribution obtained in [12, 13] in investigating the problem on a linear oscillator under 
the action of an electromagnetic field in the conditions of weak chaos. 

Let us summarize the results we have obtained above using the notions of statistical 
physics. After a lapse of time AT, that can be called the time of initial chaotization, the 



12 



investigated closed system (quantum pendulum+variable field) can be considered as a sta- 
tistical system. At that, the closed system consists of two subsystems: the classical variable 
field (6') that plays the role of a thermostat with an infinitely high temperature and the quan- 
tum mathematical pendulum (6). A weak (indirect) interaction of the subsystems produces 
passages between nondegenerate levels. After a lapse of time t > T this interaction ends 
in a statististical equilibrium between the subsystems. As a result, the quantum pendulum 
subsystem acquires the thermostat temperature, which in turn leads to the equalization of 
level populations. The equalization of populations usually called the saturation of passages 
can be interpreted as the acquisition of an infinite temperature by the quantum pendulum 
subsystem. 

IV. ENTROPY GROWTH OF THE QUANTUM PENDULUM SUBSYSTEM. 
VARIABLE FIELD ENERGY ABSORPTION 

As is known, variation or constancy of entropy can be considered as a criterion of ir- 
reversibility and reversibility of processes occurring in a closed system. In the case of ir- 
reversible processes, during which the system tends to the equilibrium state, the entropy 
increases, while in the equilibrium state it remains constant. 

Let us use this criterion to clarify the question of reversibility for our problem. As is 
known, the entropy of an arbitrary quantum system is defined by the operator of the density 
matrix p [14] 



where the brackets < . . . > denote the quantum-mechanical averaging, while the overline 
denote the averaging over a small dispersion 5uj of the passage frequency (9). In the matrix 
form the right-hand side of formula (24) is written as 



S(t) 



< p(t)lnp(t) >, 



(28) 



S(t) 



Tr{p{t)lnp{t)) 



^p ik (t)lnp ki (t). 



(29) 



■i,k 
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If initially the system is in one of degenerate states, for example, in the pure state i/) 2n , 
then C-(0) = 1, Cf(0)=0, 



Pt(t = 0) 



V 01 / 



and therefore, by (29), the entropy is S n (t = 0) = 0. 

With a lapse of time t ^> r, as passages between nondegenerate states are completed 
and the self-chaotization condition AT5u > 2n, is fulfilled, the density matrix takes form 
(22). Then the substitution of the density matrix (22) into the entropy formula (29) gives 

T 

Sn(t ~ - > r) = ln2. (30) 

Thus, on a time interval from t — to t < AT, the entropy grows 

S n {t > r) > S n (t = 0). 

This proves that on this time interval the process is irreversible. 

Using analytical methods, we have succeeded in establishing only the asymptotic value of 
entropy. To investigate a complete picture of entropy change on a time interval < t < AT, 
we use expression (15) for p+~(t). After substituting it into the entropy formula we obtain 

S n (t) = -W-(t)ln\W-(t)\ - W+{t)ln\W+{t)\ - nF n (t). (31) 

Fig. 7 shows the entropy as a function of time constructed with the aid of (18), (31) by 
numerical methods. 

To calculate the entropy value with the lapse of one period T, we substitute matrix (26) 
into the entropy formula (29) and thus obtain 

S n (t ~ T) = InA. 

The state, in which all accessible levels of the subsystem are populated with the same 
probability 1/N (Fig. 6), is the equilibrium state. The corresponding density matrix of 
dimension iV is written as 
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Pfn(t » T) 



1_ 

N 



1 . 
1 . 
1 



. 
. 
. 



(32) 



^0 I J 

After substituting p?>n from (32) into(29), we obtain the maximal entropy value on time 
intervals t >T 

S n {oo) = S n (t > T) = InN. (33) 

Thus we see that the entropy constantly grows up to value (33) and after that it stops to 
grow. 

Let us now calculate the energy mean of the quantum pendulum subsystem. It is obvious 
that for the average energy of the subsystem we can write 

N 

E(t) = E +J2Pn(t)E n , (34) 

n=l 

where E a is the initial energy value defined by the initial (t = 0) population of the levels, 
P n (t) is the probability that the n— th level will be populated at the time moment t, E n is 
the energy value in the n— th state defined from the area of nondegenerate states, In Fig. 
8 we see that the subsystem energy first grows and then becomes constant. This result 
can be explained if we take into account the time dependent trend of population changes 
which is defined by the recurrent relations (27) or Table 1. At the beginning the subsystem 
absorbs the field energy (6') and, in doing so, performs "indirect" passages between energy 
levels mostly in the upward direction. Upon reaching the equilibrium state, in which the 
subsystem is characterized by the equalization of level populations, it stops to absorb energy. 
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V. CONCLUSIONS. ANALOGY BETWEEN THE CLASSICAL AND THE 

QUANTUM CONSIDERATION 



The classical mathematical pendulum may have two oscillation modes (rotational and 
oscillatory) which on the phase plane are separated by the separatrix (see Fig.9,a). 

On the plane (E, I) the quantum mathematical pendulum has two areas of degenerate 
states - G_ and G+. Quantum states from the area G_ possesses translational symmetry 
in the pendulum phase space. These states are analogous to the classical rotational mode. 
Quantum states from the degenerate area G + possess symmetry with respect to the equilib- 
rium state of the pendulum and therefore are analogous to the classical oscillatory state [10]. 
On the plane (E,l), the area of nondegenerate states G, which lies between the areas G_ 
and G + , contains the line E = I corresponding to the classical separatrix (see Fig.9,b). If the 
classical pendulum is subjected to harmonically changing force that perturbs a trajectory 
near to the separatrix, then the perturbed trajectory acquires such a degree of complexity 
that it can be assumed to be a random one. Therefore we say that a stochastic motion layer 
(so-called stochastic layer) is formed in the neighborhood of the separatrix [7] (see Fig. 10, a). 

In the case of quantum consideration, the periodic perturbation (6) brings about passages 
between degenerate states. As a result of repeated passages, before passing to the area G the 
system gets self-chaotized, passes from the pure state to the mixed one and further evolves 
irreversibly. While it repeatedly passes through the branch points, the redistribution of 
populations by the energy spectrum takes place. Only the levels whose branch points satisfy 
condition (5), participate in the redistribution of populations (see Fig. 10, b). 
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Subscripts to Figures 

Fig.l A fragment of the parameter-dependent energy spectrum of the quantum mathe- 
matical pendulum (1). 

Fig. 2 Time-dependence of the diagonal matrix element W+(t) of the density matrix 
(15), constructed by means of formulas (15), (18) for the parameter values u = 1/t = 
1, C+(0) = 1, C~(0) = 0. As clearly seen from the figure, the higher the dispersion value of 
the parameter 5u, the sooner the stationary value W+(t > T ~ ^) = \ is achieved. 

Fig. 3 The vanishing of nondiagonal matrix elements of the density matrix (15) with a 
lapse of time t > T while the system remained in the degenerate area G_. The graph is 
constructed for the parameter values uo = 1/r = l,C+(0) = 1,C~(0) = with the aid of 
formulas (15) and (18). 

Fig. 4 A fragment of the energy spectrum depending on the slowly changing parameter 
(4) of the quantum mathematical pendulum (6). With a lapse of time t >T the stationary 
state is achieved, for which all levels are populated with an equal probability. 

Fig. 5 Results of numerical calculations performed by means of recurrent relations (27). 
Formation of statistical distribution of populations of levels P(n) with a lapse of a large 
evolution time t 1000T of the system. The result shown in this figure corresponds to the 
case for which the level population creeping is not restricted by condition (5). 

Fig. 6 Results of numerical calculations performed by means of recurrent relations (27). 
With a lapse of a large time interval t ~ 1000T the formation of stationary distribution of 
populations among levels takes place. By computer calculations it was found that in the 
stationary state all N levels satisfying condition (5) were populated with equal probability 
1/N. 

Fig. 7 The entropy growth graph constructed with the aid of expression (31), using nu- 
merical methods for the parameter values u — 1/r = 1. 

Fig. 8 Time-dependence of a mean energy value of the quantum mathematical pendulum 
subsystem, constructed by numerical methods using formula (34). As clearly seen from the 
figure, the absorption of optical pumping energy takes place prior to reaching the state of 
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statistical equilibrium. 

Fig. 9 Analogy between the classical and quantum considerations. Unperturbed motion. 

a) Classical case. Phase plane. Separatrix. 

b) Quantum case. Specific dependence of the energy spectrum on the parameter (Math- 
ieu characteristics). Degenarate G± and nondegenerate G areas of the spectrum. 

Fig. 10 Analogy between the classical and quantum considerations. Perturbed motion. 

a) Classical case. Stochastic trajectories in the neighborhood of the separatrix form the 
stochastic layer (cross-hatched area). 

b) Quantum case. The mixed state was formed as a result of population of nondegenerate 
levels situated on both sides of the classical separatrix. 
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Subscripts to Tables 

Table 1. Evolution of populations of nondegenerate levels. 

This table is a logical extrapolation of the analytical results obtained in Subsection 3. It 
shows how the population concentrated initially on one level n Q gradually spreads to other 
levels. It is assumed that the extreme upper level n + 4 and the extreme lower level n Q — 5 
are forbidden by condition (5) and do not participate in the process. 
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